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Cooperation and competition between the antiferromagnetic, d-wave superconducting and Mott- 
insulating states are explored for the two-dimensional Hubbard model including nearest and next- 
nearest-neighbor hoppings at zero temperature. Using the variational cluster approach with clusters 
of different shapes and sizes up to 10 sites, it is found that the doping-driven transition from a phase 
with microscopic coexistence of antiferromagnetism and superconductivity to a purely superconduct- 
ing phase is discontinuous for strong interaction and accompanied by phase separation. At half- filling 
the system is in an antiferromagnetic Mott-insulating state with vanishing charge compressibility. 
Upon decreasing the interaction strength U below a certain critical value of roughly U c ~ 4 (in 
units of the nearest-neighbor hopping), however, the filling-dependent magnetic transition changes 
its character and becomes continuous. Phase separation or, more carefully, the tendency towards the 
formation of inhomogeneous states disappears. This critical value is in contrast to previous studies, 
where a much larger value was obtained. Moreover, we find that the system at half-filling undergoes 
the Mott transition from an insulator to a state with a finite charge compressibility at essentially 
the same value. The weakly correlated state at half-filling exhibits superconductivity microscopi- 
cally admixed to the antiferromagnetic order. This scenario suggests a close relation between phase 
separation and the Mott-insulator physics. 

PACS numbers: 71.10.-w,71.10.Hf,74.20.-z 



I. INTRODUCTION 

One of the standard models in theoretical studies 
of high-temperature superconductors is the single-band 
Hubbard model. Substantial progress in the understand- 
ing of its ground-state properties in two dimensions has 
been achieved by applying dynamical quantum-cluster 
approaches 1 -, such as the dynamical cluster approxima- 
tion 2 -, the cellular dynamical mean-field theory 3 -^ or the 
variational cluster approach (VCA)£. Several studies pre- 
dict (i-wave superconductivity at low£ or zero tempera- 
for intermediate interaction strengths of 
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the order of the free bandwidth. 

Experimentally, it is well known that cuprate-based 
high-Tp compounds at low doping concentrations tend to 
form charge and spin inhomogeneities, such as strip o 12 ' 13 
or checkerboard modulation o 14 ' 15 : 16 ' 17 . An unbiased the- 
oretical study of inhomogeneous phases is hardly possi- 
ble within quantum-cluster approches for the presently 
accessible cluster sizes. Furthermore, for a realistic 
modelling of charge inhomogeneities additional non- 
local interaction terms should actually be taken into 
accoun t 18 ! 19 ' 20 . However, by searching for phase sepa- 
ration, an overall tendency towards the formation of in- 
homogenous phases can easily and reliably be detected 
even if only homogeneous solutions are allowed within 
the cluster mean-field calculation. 

Phase separation (PS), i.e., separation into two ho- 
mogeneous macroscopic regions with different thermody- 



namical properties, is the most simple kind of inhomo- 
geneity. It has been argue d 18 ' 20 that phase separation oc- 
curring in Hubbard- or t- J- type models with short-range 
interactions only, can transform to microscopically inho- 
mogeneous structures, such as stripes, when long-range 
repulsions were included^. For this reason, the investiga- 
tion of PS in the Hubbard model is a reasonable starting 
point for an understanding of inhomogeneous phases in 
the high-Tc cuprate materials. 

The occurrence of phase separation in the Hubbard 
model has been investigated in the past using several 
different techniques. Calculations for finite clusters us- 
ing quantum Monte-Carlo technique s] 21 ' 22 found no evi- 
dence for PS in the Hubbard model with nearest-neighbor 
hopping only. Different results have been obtained in 
infinite dimensions within dynamical mean-field theory 
(DMFT) where PS has been reported^i^. PS has 
also been discussed in the context of marginal quantum 
criticalit y 26 ! 27 , where PS occurs at the first-order side 
of the marginal quantum critical point (MQCP). For 
the two-dimensional model including hopping between 
next-nearest neighbors, calculations within the dynami- 
cal cluster approximation yield phase separation in the 
paramagnetic stated 8 -. VCA and cellular DMFT studies 
predict phase separation between a phase with long-range 
antiferromagnetic (AF) order at low doping and a super- 
conducting (SC) state at high doping 8 -^ ' 10 ' 11 . 

The new aspect of our work is to systematically inves- 
tigate the fate of the phase-separated state within VCA 
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when decreasing the interaction strength. We demon- 
strate by using as reference systems a variety of cluster 
sizes up to 10 sites, that the tendency towards PS is lost 
at small U of the order of U c ~ 4. Above of this value 
we found strong evidence from this systematic study that 
the inhomogeneous state is indeed present in the thermo- 
dynamic limit. 

In previous cellular DMFT workii, a transition which 
is in certain aspects similar to ours was found around 
U c f» 8, using 2x2 clusters only. According to previous 
studie a 29 ' 30 i 31 it is known, however, that a physical tran- 
sition for the 2D Hubbard model is taking place at much 
smaller values, U « 4. Here, the two separate energy 
scales U and J are eventually merging, as examplificd 
by the transition from two energy bands (coherent low- 
lying and incoherent Hubbard bands) to just one single 
band. A motivation for our study was partly to inves- 
tigate the relation between the metal-to-insulator tran- 
sition and PS. Most strikingly, we find that the transi- 
tion to a metallic state at half-filling occurs at essentially 
the same value of U, below which PS disappears. This 
strongly suggests a definite relation between PS and the 
Mott insulator, a relation that was previously speculated 
abouii 3 ^. 

The paper is organized as follows: In Sec. UH we intro- 
duce the model and briefly review the variational cluster 
approximation. In Sec. IIII1 our results are presented and 
discussed. The conclusions are summarized in Sec. IIVI 



II. THEORY 

Using standard notations, the Hubbard model is given 
by H = H (t) + Hi where 

H (t) = -*„.„. C L c ja ~ tn.n.n. ^ C h C j<T W 

Hi = uJ2 n n n n- (2) 



H' = H (t r ) + H x . Within the VCA one takes as a ref- 
erence system a lattice split up into isolated clusters of 
a given size. Thereby, the effects of short-range corre- 
lations on the self-energy are included on a scale given 
by the cluster extension. Trial self-energies £ = S(t') 
are varied by varying the parameters t' . Inserting the 
trial self-energy into the self-energy functional generates 
a function f2(t') = fi[53(i')] the stationary points of which 
we are interested in. It can be shown^i 3 - 3 - that 

fi(t') =n' + Trh^G^j - £(t')) _1 - TrlnGr , (3) 

where SI' is the grand potential and Gt> the Green's func- 
tion of the reference system and Go,t the non-interacting 
Green's function of the physical system. While Go,t 
is easily accessible, we calculate the reference system's 
properties using full diagonalization for small clusters, 
and the band Lanczos methods for larger ones. A de- 
scription of the numerical details can be found in Ref. [ToL 
Since we are interested in PS involving symmetry- 
broken antiferromagnetic and d-wave superconducting 
phases, our reference system includes the corresponding 
(ficticious) symmetry-breaking fields, 

H' AF = h' AF J2(nn-n i i)e i( * R < (4a) 
^sc = ^scE ! | (ClTC -'' i+h - C - ) (4b) 

ij 

where h' AF and h' sc are the strengths of a staggered mag- 
netic and of a nearest-neighbor <i-wave pairing "Weiss" 
field, rcspectiely. Furthermore, Q = (ir, it) is the AF 
wave vector, and rjij denotes the d-wave form factor 
which is equal to +1 (—1) for nearest-neighbor sites with 
Ri — Rj in x (y) direction. In addition, an on-site po- 
tential is included in the set of variational parameters to 
ensure a thcrmodynamically consistent determination of 
the average particle number—. Note that all three varia- 
tional parameters couple to one-particle operators only. 



The operator c\J creates (annihilates) an electron with 
spin a at the site i, and rii a is the corresponding oc- 
cupation number operator. We consider both hopping 
t n .n. along nearest-neighbor bonds (ij) as well as hopping 
tn.n.n. along next-nearest-neighbor bonds ((ij}). U is the 
local Coulomb interaction. We set the unit of energy by 
t n .n. and choose t n . n .n. = — 0.3i n . n . throughout the paper 
which is a realistic value for the cuprate materials. 

The main idea of the variational-clustcr approximation 
(VCA)^!^ 3 - is to consider a "reference system" to span a 
space of trial self-energies among which the self-energy 
that describes best the physics of the infinite-size lat- 
tice model is obtained via a dynamical variational prin- 
ciple <5f2[S] = 0. Here fl stands for the grand poten- 
tial. The reference system is given by a Hamiltonian H' 
with the same interaction part Hi as the physical sys- 
tem but with modified one-particle parameters i', i.e., 



III. RESULTS 

For the calculations we concentrate on the hole-doped 
side of the phase diagram since it has been shown that 
there PS is much more pronounced than in the electron- 
doped case&Siiii. The occurence of PS can be best in- 
ferred from the dependence of the chemical potential on 
the particle number. For a physical system in thcrmody- 
namical equilibrium the charge susceptibility k — dn/dfi 
must be non-negative. Hence, k < indicates a thcrmo- 
dynamically unstable phase. At a fixed average density 
n lying in the instability region n\ < n < ri2, the free 
energy can be reduced if the system develops two spa- 
tially separated homogeneous phases, one with a fraction 
x = ™^Zn ± °f particles at the density n\ < n and an- 
other one with the fraction 1 — x at ri2 > n, rather than 
having a single homogeneous phase. The boundaries ri\ 
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FIG. 1: (Color online) Chemical potential fi as function of the 
average density n, calculated with a 2 x 2 cluster as reference 
system for different values of U. From top to bottom: U — 12, 
8, and 4, respectively. n\ and rii denote the boundaries of the 
instability (PS) region. Black solid lines: Pure SC phase. Red 
dashed lines: AF+SC mixed phase. Dotted horizontal line: 
critical chemical potential (where applicable). 



and ri2 of the instability region can be obtained by a 
Maxwell construction (see Fig. [J) , in close analogy to a 
gas-liquid system. After the Maxwell construction, the 
physical chemical potential is independent of n between 
ni and ni. 

In Fig.Q] we show results for the interrelation of the 
chemical potential [i and the average density n as a func- 
tion of Coulomb interaction U. The calculations have 
been done with a 2 x 2 cluster as a reference system. In 
the strong coupling regime, U = 12, we find that the 
chemical potential [i as function of the density n shows 
a nonmonotonic behavior (and thus k < 0) as mentioned 
above. We conclude that for this coupling and within 
the precision set by the 2x2 cluster reference system, 
the results imply PS into an AF+SC mixed phased at 
low doping and a purely SC phase at higher doping. The 
critical chemical potential where the two phases coexist, 
is indicated by a horizontal dotted line in Fig.[TJ This 
behavior is very similar to the previously reported one 
for U = 8*^. For comparison, results for U = 8 - as 
published in Ref. - are shown in Fig.Q] in the middle 
panel. 

The picture changes even qualitatively when going to 
weak interactions, e.g., U = 4. As shown in the lower 
panel of Fig.Q] the chemical potential as function of n 
shows a positive slope for all dopings, without an inter- 
mediate region with negative sign. We conclude that for 
weak interaction there is no tendency to the formation of 
inhomogeneitics, even for the smallest cluster we used for 




FIG. 2: Clusters used as reference systems in this study. The 
"a/TO x VW" and "VS x y/E" clusters are marked explicitly 
in the figure. 
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FIG. 3: (Color online) Characteristic energy A^i = fj,* — fi c 
(top) and critical densities n\ and ri2 (bottom) as functions 
of U, calculated using the clusters shown in Fig.f5]as reference 
systems, 2x2 (black, circles), \/8 x y/8 (blue, squares), 4x2 
(red, stars), 5x2 (magenta, triangles down), and \/l0 x \/l0 
(green, triangles up). Lower plot: Solid lines correspond to 
the critical density ni , dashed lines to the critical density ri2 . 



our calculations. Instead, the AF+SC solution is stable 
up to larger dopings, and the staggered magnetization 
vanishes continuously in the stable and homogeneous so- 
lution. 

Since the occurrence of PS is subject to rather strong 
finite-size effects^, and in order to further elucidate the 
difference between the coupling regimes, we have recalcu- 
lated the phase diagram using larger clusters with L c = 8 
and L c = 10 lattice sites as reference systems (see Fig.f2|). 
For both L c = 8 and 10 we have considered two differ- 
ent cluster geometries to estimate dependencies on the 
cluster shape. 

In addition to the boundaries of the instability region, 
ni and n, 2 , we define a characteristic energy A/i = fi*—fi c . 
Here, fi* is defined as the point where the slope of 
n(n) changes sign, and \i c is the critical chemical potcn- 
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FIG. 4: (Color online) Ground-state energy per site Eo/L 
as function of electron density n, for U = 4 (top), f7 = 8 
(middle), and U = 12 (bottom). Results are shown for the 
2x2 (black solid), 4x2 (red dashed-dotted), and 5x2 (green 
dashed) reference systems. 

tial, i.e., the point at which ^af+sc (m) crosses f2sc (/•*)• 
Equivalcntly, fi c is fixed by the Maxwell construction, 
see Fig.[TJ In the limit of infinite cluster size, the char- 
acteristic energy A/i must vanish as the reference sys- 
tem (and thereby also the original system) is solved ex- 
actly within the VCA: For L c — > 00 (and for densities in 
the exact instability region) the reference system spon- 
taneouly generates the phase-separated state, and n(n) 
becomes flat between n\ and 712, as discussed above. On 
the other hand, the difference between the critical densi- 
ties ri2 — n\ must converge to a nonzero value for L c — > 00 
whenever the system in the thermodynamic limit shows 
phase separation. Notice that the difference between the 
grand potentials f^AF+sc and Qsc becomes smaller and 
smaller with increasing cluster size. Thus, together with 
the increasing numerical effort, this makes the identifi- 
cation of A/i harder for larger clusters. Contrary, the 
slope of Q(/i), i. e. the particle density, is not affected by 
a systematic shift and, therefore, its calculation is quite 
reliable also for larger L c . 

Results for these quantities are shown in Fig. [3] The 
upper panel shows A/i as a function of the Hubbard in- 
teraction U. As expected from the discussion above, A/i 
decreases for increasing cluster size. Quite surprising, 
however, is the dependence of A/i on the interaction: 

As A/i 7^ is a finite-size effect, one could expect A/i to 
be smaller for stronger interactions since a cluster mean- 
field approach quite generally is expected to be more re- 
liable in the strong-coupling than in the weak-coupling 
regime. At least, it is known from previous calculations^ 
that for weaker interactions finite-size effects are larger 
as could be seen from the stronger dependence of the op- 
timal variational parameters on the size of the reference 
system cluster. In fact, we find for large U (see Fig. [3]) 
that Afi increases with decreasing interaction. Below 
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FIG. 5: (Color online) Order parameters as a function of 
density for U = 4 (left), [7 = 8 (middle), and U = 12 (right). 
Top: AF order parameter. Bottom: d-wave SC order param- 
eter. Results are shown for the 2x2 (black solid), 4x2 (red 
dashed-dotted), and 5x2 (green dashed) reference systems. 



U as 10, however, A/i starts to decrease and eventually 
even vanishes at some critical value U c (which exhibits 
a finite but weak dependence on the cluster size). This 
decrease indicates a qualitative change of the phase tran- 
sition with decreasing U. 

From the analysis of the critical densities we see that 
the vanishing of A/i corresponds to the disappearance 
of the phase-separated state: The critical densities ri\ 
and 712 of the two coexisting phases are shown in the 
lower panel of Fig.[3j The dependence of their difference 
712 — Hi on U closely resembles the dependence of A/x. In 
particular, for U below some U c , n\ and ri2 collapse to a 
single point, i.e., phase separation disappears. 

It is encouraging to see that, except for the vTO x vTO 
cluster, the critical densities shown in this plot depend 
only very weakly on the cluster size. The VTO x VlO 
cluster seems to behave differently. Probably, its shape 
(see Fig. ^ makes this cluster unsuitable for finite-size 
scaling. This is also confirmed by the fact that A/j, does 
not display the correct scaling behavior for this cluster. 
In addition, the 5x2 cluster with the same number of 
sites L c = 10 does display the correct L c dependence of 
A/i. Excluding the results from the vTfJ x V40 cluster, 
ni and n-i show a very weak dependence on L c . Close to 
U c , the L c dependence is somewhat stronger again which 
is not surprising, of course. 

In view of the results for different interaction strengths 
and cluster sizes we conclude that PS persists in the ther- 
modynamic limit down to a critical interaction U c . Ad- 
mittedly, there is a rather large uncertainty in the de- 
termination of U c : Estimates range between U c ~ 4 and 
U c ~ 6. Interestingly, results from the VTO x VlO cluster 
do not only show a lower critical interaction but also an 
upper critical value [/ c2 ~ 10 above which PS disappears. 
As discussed above, however, the question is whether re- 
sults for this cluster shape are reliable or not. 
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FIG. 6: (Color online) Closing of the Mott gap as function 
of U. From top to bottom: L c = 2 x 2, L c = 4 x 2, and 
L c = vlO x \/W site clusters. From left to right: U — 3 
(green), 3.5 (blue), 4 (black), and 5 (red). 



From the grand potential we can extract the ground- 
state energy of the the system by Eq = fi + /J.N. Re- 
sults for the ground-state energy per site, i. e. Eq/L, arc 
shown in Fig.[3]for U = 4, 8, and 12, resp. Although some 
(small) finite-size effects are visible, the ground-state en- 
ergy seems to be well converged in our approach. As a 
check for consistency we found that the general relation 
d(E(j/L)/dn = /i is fulfilled within numerical accuracy. 
This relation is a consequence of the Legendre transfor- 
mation Eq = fi + 1-iN and must be true in all calculations 
independent of the cluster used as a reference system. 

The filling dependence of the order parameters is dis- 
played in Fig. [5] Although some finite-size effects are vis- 
ible, the results strongly indicate that superconductivity 
persists in the thermodynamic limit. We can also ex- 
tract an interesting trend when comparing the results for 
U — 8 with those for U = 4: While the staggered magne- 
tization decreases the SC order parameter increases with 
decreasing U. At the same time the doping region where 
a magnetic solution exists, extends to somewhat larger 
dopings for smaller U. Another difference is the small 
but finite value of the dSC order parameter at half-filling 
(n = 1) for U = 4, as can be seen in the lower right plot 
of Fig. [5] This is directly related to the closing of the 
Mott-Hubbard gap, as discussed below. 

There are also some differences between U = 8 and 
U = 12, although they are considerably smaller. Tak- 
ing for instance the boundary of the AF phase, i. e. the 
density where the AF order parameter vanishes, it is at 
n ps 0.82 for U = 4, n ps 0.87 for U = 8, and n ps 0.89 
for U — 12. Also the changes of the SC order parameter 
are very small. We want to stress again that the major 
difference, however, is the absence of PS for [7 = 4, see 
Figs. [I] and gj 



To understand the occurrence of a critical interaction 
for PS as well as the finite dSC order parameter at half- 
filling, we study the n(/z) behavior close to half-filling 
as a function of U. Results are shown in Fig. [5] in the 
vicinity of n = 1. For U = 5 there is a flat n(/x) de- 
pendence for a finite range of chemical potentials, i. e., a 
vanishing charge susceptibility k = 0. This indicates the 
presence of an (AF) Mott insulating state at half-filling. 
The Mott gap shrinks with decreasing U, and eventu- 
ally the system shows metallic behavior at n = 1. This 
agrees well with previous result a 30 ' 31 ' 36 i 37 . For the small- 
est system the transition to the metal takes place at an 
interaction strength somewhat above {7 = 4, whereas for 
larger clusters the critical U is shifted to slightly smaller 
values. For the 5x2 cluster (not shown) it is slightly be- 
low U — 4, and for the \/T0 x ylO cluster the transition 
takes place slightly above C7 ~ 3.5. Notice that the value 
of U where the gap closes can also be read off directly 
from the dSC order parameter. Below the critical U, the 
system is metallic even at n = 1, and hence can also be 
superconducting at half-filling. This is seen in our cal- 
culations which give a finite order parameter at n = 1 
(see Fig. [5]) . Wc like to stress that close to half-filling the 
solution with lowest energy is always a mixed AF+SC 
one^i. We checked this by comparing this solution to the 
pure AF and SC ones. 

It is obvious that our determination of U c is not as 
precise as it can be done by other methods, e. g. the 
path-integral renormalisation group^ 7 -. Moreover, with 
the cluster sizes available in the present form of the VCA, 
it is not possible to detect whether or not we have a 
MQC P 26 i 27 with a diverging charge susceptibility k at 
the metal-insulator boundary. Nevertheless, our calcula- 
tion allows to make a qualitative connection of the MIT 
and PS. It indicates that the collapse of the upper and 
lower dopings n\ and ri2 (i. e. the disappearance of PS) 
and the closing of the Mott gap occur at almost the same 
critical interaction strength. It appears that it is impor- 
tant to have a Mott insulator at half-filling in order to 
get phase separation away from half-filling. Due to the 
limited cluster sizes available and the according limited 
accurary in the determination of critical points, however, 
this statement is somewhat speculative. Note that a rela- 
tion between phase separation and the doping of a Mott 
insulator has been discussed in Rcf. 32. There, it has 
been argued that a Mott insulator with broken symme- 
try, which is the SU (2) symmetry in our case, enhances 
the possibility of PS. 

In order to analyze the relation with single-particle ex- 
citations, we plot in Fig. [7] the Fermi surface for weak 
coupling, U = 3.5, calculated with the v8 x V8 clus- 
ter. The closing of the Mott gap is clearly manifested 
by the occurrence of both hole and electron pockets near 
(§,§) and (7r,0), respectively, in direct contrast to the 
Fermi surface for larger coupling and U = 8 where one 
has either hole or electron pockets at low doping^. The 
simultaneous occurrence of hole and electron pockets has 
also been reported recently in weak-coupling calculations 
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FIG. 7: Representation of the Fermi surface for an interaction 
U = 3.5 < U c and n = 1. Calculation using the y/8 x y/8 
cluster. Dark regions denote large spectral weight integrated 
over a small frequency window of Aui = ±0.05 around uj = 0. 
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FIG. 8: Spectral function for U = 8 (top) and U = 4 (bot- 
tom). Calculation using a \/8 x v8 cluster. The horizontal 
dashed line marks the chemical potential. 



close to half-filling^. 

Another way to see the closing of the Mott gap is to 
look at the spectral function directly, which we plot in 
Fig. [8] for U = 8 and U = 4, resp. Calculations have 
been done using the same V8 x V8 cluster as for Fig.0 
but at finite hole-doping just at the critical doping n-i in 
the pure SC phase. From this figure, it is obvious that 
for [7 = 8 one has a quasi-particle band at the Fermi 
level, well separated from the upper Hubbard band. For 
U = 4, however, the upper and lower Hubbard bands 
merge, and closely resemble the band structure of the 
free system, except for the superconducting gap around 
(tt.0). 

Related ideas of a change of the type of the phase tran- 
sition as function of the interaction strength have been 



reported quite recently based on cellular DMFT calcu- 
lations by Capone et aliii for the two-dimensional Hub- 
bard model with t n . n . n . = 0. Capone et al. inferred a 
critical interaction of approximately U c ~ 8 below which 
the doping-dependent transition to a paramagnet is con- 
tinuous and above which there is a first-order transition 
accompanied by phase separation. The U c reported for 
tn.n.n. = and using cellular DMFT with 2x2 clusters 
is somewhat larger as in our study. A qualitative differ- 
ence, however, is that we always find a mixed AF+SC 
state at low doping whereas the authors of Ref. [ll| find 
a pure AF phase for strong interaction close to half- 
filling. One should note, however, that finding the sta- 
ble solution with lowest energy can be a difficult task 
within quantum-cluster approaches. In Ref. [0 the cal- 
culation starts from pure AF and SC solutions and a 
mixed AF+SC solution is searched for by applying small 
perturbations to the pure ones. This procedure does not 
necessarily lead from one local solution to another, espe- 
cially if the two are well separated in parameter space. 

Finally, we would like to comment on the relation of 
our results to the t-J-modcl. Translating our results of 
Fig. [3] to J oc jj, we argue that PS for small J (large 
U) should be weak, or eventually absent in the thermo- 
dynamic and U — > oo limit, since Afi decreases signifi- 
cantly with increasing U. However, from our results it is 
not possible to deduce a definite critical J c below which 
PS should be absent. Nevertheless, the spectral function 
and in particular the low-lying quasi-particle band in our 
calculations is in well agreement with the low-lying bands 
in the t- J- models 



IV. CONCLUSIONS 

To summarize, we have studied the occurrence of phase 
separation in the two-dimensional Hubbard model at zero 
temperature depending on the strength of the Hubbard 
interaction U. We have employed the variational-cluster 
approach using clusters of different shapes and sizes up to 
10 sites. Our results show that the nature of the doping- 
dependent transition from the antiferromagnetic and su- 
perconducting state to a non-magnetic and purely super- 
conducting state changes from discontinuous to continu- 
ous when going from the strong- and intermediate- to the 
weak-interaction regime. Below a critical value for the 
interaction strength no phase separation can be found. 
Above the critcal value, there is a clear discontinuity in 
the density as a function of the chemical potential with 
a small finite-size error only. This is a strong indication 
that PS, or at least a phase with mesoscopic inhomo- 
gencitics, should persist in the thermodynamic limit, as 
long as the interaction is sufficiently strong. 

We also studied the interaction-driven Mott transi- 
tion at half-filling. For weak interactions (and t n . n . n . = 
— 0.3f„.„.) the system is in a mixed AF+dSC state with 
a finite charge susceptibility. The transition to the AF 
Mott insulator with vanishing charge susceptibility takes 
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place at U ~ 3.5 — 4. As the Mott transition at half- 
filling takes place at the same interation strength where 
phase separation appears in the doped case, one might 
speculate that these two phenomena are closely related. 
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